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Abstract 

We study neutrino oscillations in space within a realistic model in which both 
the source and the target are considered to be stationary having Gaussian- 
form localizations. The model admits an exact analytic solution in field theory 
which may be expressed in terms of complementary error functions, thereby 
allowing for a quantitative discussion of quantum-mechanical (coherent) ver- 
sus statistical (incoherent) uncertainties. The solvable model provides an 
insightful framework in addressing questions related to propagation and os- 
cillation of neutrinos that may not be attainable by the existing approaches. 
We find a novel form of plane-wave behaviour of neutrino oscillations if the 
localization spread of the source and target states due to quantum mechanics 
is of macroscopic size but much smaller than neutrinos' oscillation length. 
Finally, we discuss the limits on the coherence length of neutrino oscillations 
and find that they mainly arise from uncertainties of statistical origin. 
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It was thirty years ago when Pontecorvo suggested |] that if the known light neutrinos 
u e , and v T have non-zero masses and mixings, they may oscillate to each other in the 
very much the same way as the K° and K° mesons do 0. This mechanism has also been 
invoked to explain the energy deficit between the neutrinos produced in the sun or earth's 
atmosphere and those detected on earth's ground || . In particular, according to the recent 
results of the Super-Kamiokande collaboration M, the experimental data may comfortably 
be explained by naively assuming u^-to-v? oscillations, when compared to the non-oscillation 
scenario of the Standard Model. For neutrino energies E much larger than their respective 
masses m 8 (with i = 1,2,3), the classical formula |T[ governing the probability of v a -to-Up 
oscillations as a function of the distance / from the source reads 
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W) = E iVaifm 2 + 2£ XetYaiV&V^Vpu) cos( lm ^ ~ mfcJ ) , (1) 
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where V a j is a leptonic mixing matrix analogous to the Cabbibo-Kobayashi-Maskawa mixing 
matrix of the quark sector. For simplicity, we have also neglected possible effects of CP 
violation in Eq. (Jl]). Up to now, many authors |5|-|iTf have discussed the conditions under 
which the above naive cosine formula of oscillations between neutrinos or other particles 
may or may not be valid. In a large number of papers fl|6|,|^ Jl^ , [l3|JT7| , the derivation of Eq. 
(0) has been based on using neutrino wave packets. Within other more field-theoretically 
oriented approaches |7|, pl] , p!2| , p^| , neutrinos are treated as virtual intermediate states which 
are produced by some initial states and are only observed through their interaction with 
the target state in the detector. The latter considerations may also be related to those of 
the old-fashioned S-matrix approach |2,18]. However, in all different treatments present in 



the literature, many approximations were necessary to arrive at the known formula ([!]) of 
oscillations. 

In this paper, we shall study particle oscillations in space within a model in which both 
the source and target states are taken to be stationary having Gaussian-type broadenings. 
For our illustrations, we first consider the oscillation of two scalar particles whose dynamics 
may well be determined by the relevant propagators at the second order of electroweak 
interactions. The discussion is then extended to the case of fermions, i.e., that of neutrinos. 
Since we assume a Gaussian-form localization for the production and detection interactions, 
the model admits an exact analytic solution within the framework of field theory and hence 
different sorts of approximations can directly be controlled for their validity. Thus, the 
solvable model offers an important insight into the profound issue of particle oscillations, 



thereby complementing related approaches discussed in [0,|TT].|12],|14[ . In particular, we find 
that the oscillating pattern of neutrinos depends crucially on the coherent broadening of 
production and detection points. As long as the coherent broadening effects lie within the 
microscopic range, one has the usual picture of particle propagation through spherical waves 
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in the three-dimensional space, as the particles are emitted from a well-localized source, e.g., 
neutrinos from the sun. If such a broadening, however, happens to be of mesoscopic size but 
still is much smaller than the oscillation length of particles, then their propagation proceeds 
via plane waves and their oscillation length depends explicitly on the three momentum 
direction of the initial states at the source. 

Let us consider the 2^2 scattering process with two scalar particles as intermediate 
states which may oscillate to one another. The production and detection of the intermediate 
scalars occurs through some asymptotic states at the space points x and y, respectively. 
For thermal-source situations under study |IB| , the time elapsed between production and 



detection cannot be measured directly. Therefore, it is legitimate to assume that to a rather 
good approximation, the interactions taking place in the source and target are stationary. 
The two points x and y are then macroscopically separated by a distance I = \x — y\. 
For example, for solar neutrinos I is the distance between sun and earth, for atmospheric 
neutrinos I is the distance between earth's atmosphere and detector at the ground, etc. 
Furthermore, we take the spatial localizations of the interactions at the source and target 
to be Gaussian functions peaked at x and y with dispersions of cubic form, i.e., 5x % = 5x 
and by 1 = 5y. 

With the above considerations, the amplitude for a — > (3 transition is given by 

Tp a {k,y,5y;p,x,5x) = Af [ d 4 xd 4 y e -(---)7^ 2 e -(i?-s)W e -ipx+iky 
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where p = P1+P2 and k = k\ + k2 are the total momenta of the initial and final states, respec- 
tively. The overall proportionality factor M contains coupling constants and other irrelevant 
multiplicative terms; it drops out after the probability rate is normalized. Furthermore, 

Of = O aj O Pj (3) 

is the usual combination of elements of the mixing matrix O. The 2x2 orthogonal matrix 
O diagonalizes the mass matrix of the two-scalar system, and relates the weak to the mass 
eigenstates. Obviously, Y^jOj = S a p. Here, we should also remark that Lorentz invariance 
of the amplitude @ is not manifest; it is broken by the space-dependent Gaussian functions. 
However, one can always recover the manifest form of Lorentz invariance by considering 
instead the Lorentz-invariant exponent Sk^(x — x) M for the Gaussian function with 6k^ = 
(1/Sx° ,1/Sx l ) and likewise for the y coordinate. The stationarity condition implies that 
5x° 3> 5x l and Sy° 3> Sy l . For our inertial frame of particle production and detection, we 
take 5x°, 5y° — > oo. 

Since the interactions at x and y do not display any time dependence in the stationary 
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limit discussed above, the integration over the times x° and y° can now easily be performed. 
The transition amplitude Tp a then takes the form 

-° — rP") ( rl 3 rrl 3 il p -{x-x) 2 /Sx 2 „-{y-y) 2 /Sy 2 ipx-iky 
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where q 2 = E 2 — m 2 . The result of the time integrations is that energy conservation holds 
strictly both at production and detection vertices, i.e., k° = p° = E, as can readily be seen 
from the delta function on the RHS of Eq. (|4]). 

It is now important to notice that under spatial displacements, e.g., x — > x + a and 
y — > y + a, Tp a (x,y) transforms into e m ^ k ~^Tf3 a (x,y). In other words, the amplitude is 
frame-independent up to an unobservable phase a(k — p). Employing this fact, it proves 
convenient to redefine the amplitude into a manifestly frame-independent form 

Tp a (k,p,l,5y,8x) -> e - l P x+lk y Tp a (k,y,5y;p,x,5x) 

= Af5(k° -p°) J d 3 xd 3 ye-* 2/Sx2 e"^ 2 e^-^ 

2 a/3 f d 3 q Q-^-y- 1 ) 
x £ °1 I (2^8 q i-\q\2 + i£ ' (5) 



with I = y — x. It is now not difficult to carry out the Gaussian integrals over x and y (see 
also Eq. (|A.1| ) in Appendix A). In this way, we obtain 

2 a r d 3 f7 pSx 2 (p-q) 2 /4 5y 2 (k-g) 2 /A iql 

T, a (k,p, 1, 6y, 5x) = AfS(k° - p°) £ Of J ^ <$ - \tf + ie ' (6) 

with M = 7i 3 5x 3 5y 3 N. In Appendix A we show that the very last three-dimensional integral 
over q in Eq. @ can be solved exactly within a class of error functions. Taking this very 
last fact into consideration, we arrive at the following analytic expression: 

T f3a (k,p,l,5y,6x) = -M5(k°-p ) — £ e -\^ 2 M 2 +q 2 )-\SyH\k?+^) 

87r|L| j = i 

— * — * 

x >^ Erfc( - '-Slq, - H) _ e -**ffl Erfc ( _ »^ + H) ] , (7) 

where Erfc(z) is the complementary error function defined in the appendix, 8l 2 = 5x 2 + Sy 2 
and 

L = I - l -5x 2 p- l -5y 2 k. (8) 

The norm of the complex vector L is understood to act only on the vectorial space, i.e., 
\L\ = \Jl~ 2 . Equation (|7]) is the major result of this paper, from which any behaviour of 
neutrino oscillations in space for different kinematic conditions may directly be inferred. 
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On physical grounds, one generally expects that k — p = 0(l/5x,l/5y). Therefore, 
without sacrificing any of the physical features of particle oscillations we wish to study, 
we can explicitly impose a kind of three-momentum conservation for the initial and final 
states of the scattering in Eq. (R), i.e., k = p This leaves the three- momentum of 



the oscillating particles still unspecified. In fact, the very same result would have been 
obtained, if we had started with a Gaussian interaction of the form exp[— (x — y — l) 2 /5l 2 ], 
with 5l 2 = 5x 2 + 5y 2 and integrated out first the frame-dependent components (x + y)/2. 
However, we should stress that our conclusions would not change even if we considered 
the complete but more lengthy expression (0). Under these assumptions, the probability 
amplitude simplifies to 

% a (k,l,5l) = -M5^\k-p) ^ =r V Of e -3 5/2 ^l 2+ ^ \ e ^\H^\ 

8tt \l - %5l 2 h\ J^x L 

xErf C ( - i« B - £=M) _ ^ir-i^-i Erfc( _ y lq . + lizMi) ] . (9) 

Notice that Af must be redefined accordingly for reasons of dimensionality. We can now 
determine the probability for the state a produced at a distant well-localized area, e.g., the 
sun, to observe the state (3 at a macroscopic distance /. The normalized probability is given 
by 

d p (k i m \Tp a {k,l,5l)\ 2 

where T denotes the phase space of the final states. In Eq. (|10|), we have not yet included 
uncertainties of statistical origin which can be added incoherently (we defer the discussion 
to the end of the paper). In the following, we shall focus our attention on coherent quantum- 
mechanical uncertainties. 

It is very instructive to discuss the following two limiting cases: (i) 5l 2 \k\ I and (ii) 
5l 2 \k\ ^> I. Let us first consider case (i). We can expand the norm of the complex vector in 
Eq. (H) in terms oi 5l 2 \k\, i.e. 

ir-^jEi = ,_•«» f +0 («J!!). (id 

To leading order, the transition amplitude then takes on the known oscillatory form 

T Pa {k,T,8l) = -M8^{k-p) -L V Of e ~W-*W e *oi + o(^P~) . (12) 

Ant ~^ J v r ' 

The fc-dependent exponential factor in Eq. ( |T^) controls the flow of the three-momentum 
between the asymptotic states and the intermediate particles, i.e., that of neutrinos. In 
fact, the three-momentum is conserved in the production and detection vertices up to an 
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error of order 1/51. In the limit 51 — *> 0, any information about the three-momentum of the 
oscillating system is completely lost, e.g., the total three- momentum \k\ of the detecting 
particles can take any possible value. However, in experiments the distribution over k can 
in principle be measured, so physically 51 may be small but not zero. How small or how big 
51 could be is a puzzling issue that will be discussed below. 

After integrating over the three-momentum phase space V in Eq. (PH|), we find the well- 
known formula of particle oscillations 

P a ^(E,l,5l) = £ Ofe^ 1 ] + o(5l( qi -q 2 ),°-^) . (13) 

The above formula is valid under the following two conditions: 

2ir 

I. 51 < L osc = , (14) 

qi ~ <?2 

II. 51% < I . (15) 

The first requirement reflects the fact that the coherent uncertainty in the position 51 for the 
source and target states must be much smaller than the oscillation length L osc in order to have 
a non-vanishing oscillating pattern. Within the field-theoretic framework, this constraint 
has first been discussed in Ref. [0 and further elaborated on in |TTlJl^ , [T4]JT5l . Note that a 
constraint analogous to Eq. (|T4| ) may also be obtained in the wave-packet treatment if 51 
is to be interpreted as the size of the neutrino wave packets P,|6], [10| , |13| ) |I^ , [17| . Nevertheless, 
the new requirement in Eq. (|i"5|) also affects the form of particle oscillations and deserves 
high attention as well. In particular, it is interesting to notice that one can satisfy condition 
I by grossly violating condition II. To give an example, suppose that 51 = 1 cm, Am 2 = 
m\ — m 2 , = 10~ 4 eV 2 , and the energy E = 1 GeV (1 GeV ~ 5 x 10 13 cm -1 ) which is typical 
for atmospheric neutrinos. Then, one has L osc « E/Am 2 « 10 8 cm, certainly much larger 
than 51. On the other hand, one finds that 5l 2 q\ ~ (1 cm) 2 x 1 GeV k 5 x 10 13 cm is 
comparable to the distance between sun and earth! 

To elucidate further the implications of the new condition II for particle oscillations, we 
shall now discuss the behaviour of the probability amplitude for the limiting case (ii) when 
5l 2 \k\ ^> I. In this limit, the norm of the complex vector I — ^5l 2 k may be expanded as 
follows: 

\l--5l 2 k\ = - -5l 2 \k\ + ^ + Of— . (16) 
' 2 1 2 1 1 \k\ ^5l±\k\ 2) y 1 

Substituting the last expression into Eq. (|9|) and employing the asymptotic expansion of 
Erfc(z) for large values of z in Eq. ( |A.15|) yields 

T 0a {k,l,5l) = -M6W(k-p) V Of e-^ 2 d fe >*) 2 e<**T/l*l + o( ^ 

2n5l 2 \k\ h 3 K SlMk\ 



(17) 
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Because of the absence of the l/l dependence of Tp a in Eq. (|17])> the propagation of the 
intermediate states turns over to plane waves depending explicitly on the three momentum 
k. This intriguing plane-wave behaviour of the amplitude may be attributed to the fact that 
non-localization effects of the source and target quantum states within a finite space volume 
fl9fl get coherently amplified for large (mesoscopic) values of 51. Geometrically, one may 
attempt to visualize this counter-intuitive result of the behaviour of particle propagation as 
follows. It appears that the intermediate states are emitted and/or detected by means of an 
'antenna' of macroscopic size Sl 2 \k\ instead of 51 that one would have naively expected. So, 
if an observer were close to such an 'antenna' emitting neutrinos, say, in a distance much 
smaller to its size, she/he would find in principle that these particles are emitted by plane 
waves, whereas she/he would only recover the usual form of propagation through spherical 
waves if she/he were at a distance much bigger than the size of the 'antenna.' Nevertheless, 
we still feel that a deeper understanding of this quantum-mechanical phenomenon would be 
very useful. 

One might now worry that averaging effects due to the phase-space integration of the 
final states would cancel the factor proportional to kl which appears in the oscillating part of 
the probability in Eq. (|i~7|) . However, this is not quite true, since the square of the transition 
amplitude Tp a (k, I, 51) does not depend on the three-momentum difference k_ = ki — k 2 , and 
hence the phase-space integration over T has the effect to replace simply the oscillatory 
factor kl with pi. To make this explicit, we note that one can always make the following 
variable substitutions in the phase-space integral: 

d 4 k x d*k 2 5 + (kl-m\)5 + (kl-ml)5 {i \k-p) \T Pa {kJ,5l)\ 2 

= ~ J d 4 k d 4 k. 5 + [ {k + 4 fc ~ )2 - m\] 5 + \ {k ~ A k - )2 - ml] * (4) (fc ~ p) \T 0a (k, 1 5l)f 
1 



^\T Pa (p,l,5l)\ 2 I d 4 k- 5 + [ (P + A k ~ )2 - m 2 ] 5+ [ (P - ml] , (18) 
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where 5 + {kf — mf) = 9(k®)5(k 2 — m 2 ) (i = 1, 2). The integral in the last step of Eq. ( |I~8"D 
is just an overall normalization constant and cancels in Eq. (|i~0D. Taking this into account, 
the probability of particle oscillations reads 

P a ^(E,p,l,6l) = |E Of e^^\ 2 + o(5l( qi -q 2 ),^-) . (19) 



From the above formula, it is obvious that the oscillation length as seen by the detector is 
strongly correlated with the three-momentum p of the initial states producing the mixed 
particles. In principle, the two different production mechanisms as well as their predictions 
obtained by Eqs. ( |T3| ) and ([19]) could be distinguished if we were able to measure the spatial 
dependence of the oscillation length by moving the detector around the source. In addition, 
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we must assume that one somehow knows the initial momentum p of the production mecha- 
nism and controls well both coherent (5y) and incoherent uncertainties at the detector. Even 
though this seems to be a rather formidable task, the actual size of 51 dictates the nature 
of particle propagation within the solvable model. So far, many qualitative estimates exist 
for 51. However, there has not been yet a rigorous method to evaluate this quantity from 
first principles. It is obvious that 51 depends decisively on the details of the experiment. 
According to ||,|13|,[1(],|16|] , the most natural choice for the scale of 51 for neutrinos coming 
from the sun lies between the size of nucleus, i.e. several fm's, and 10~ 7 cm. The latter 
originates from effects due to thermal collisions of electrons which results in a reduction of 
the wave-packet size of the production system. Therefore, it is most likely to expect that 
solar neutrinos have a spherical-wave propagation. 

We shall now turn to the case of oscillations of two fermions such as neutrinos. Our 
starting point is a transition amplitude analogous to Eq. (|) 

Tp a (k,pX5y,5x) =N 5(k° - p°) Jf{h)up{k 2 )^P L J d 3 xd 3 ye 

d 3 q (>y E - 7g + rrij) e -^-£*- r ) 



-x 2 /Sx 2 Q—y^/Sy 2 ^ipx—iky 



ual 3 f d^q ^E- iq + m 3 e~^-^> 
§^ ij^f d-W + ie 7^«a(P2)4 (Pi 



X 

= Af5(k°-p°) r f (k 1 )up(k 2 ) 1 ^P 1 



X 

3=1 



2_ r d 3 q e - fa2(?7 "^ 2/4 e - S y 2 (k~q) 2 /4 e iql 

(2tt) 3 q 2 -\q\ 2 + ie 



E V? (J°E - 7 V f ) / 



x wMjrfa) . (20) 



where p = p\ + p 2 , k = k\ + k 2 , Pl = (1 — 7s)/2, E = k° = p° (qj = ^E 2 — mj) and 

vf = v a3 v; 3 . (21) 

Specifically, Eq. ( p0|) describes the amplitude probability for producing the neutrino flavour 
a at point x through the standard V — A weak current interaction J^(pi) and detecting the 
neutrino flavour (5 through a similar V — A weak current interaction Jf(k\). The neutrino 
flavour states a and j3 may be identified by the charged leptons l a and lp which accompany 
the neutrinos in the production and detection vertices, respectively. 

The analytic form of the transition amplitude describing neutrino oscillations is given by 

T Pa (k,pX5y,5x) = - M 5(k° - p°) — £ e^ 5 ^**®-*^®'*® 



8tt L\ 



3=1 



x ( J l }{k 1 )u p {k 2 ) 1 ^P L E 1 ° lv u a {p 2 )J^{p x ) (e^^Erfcz^ - e'^Ericz^ 



- J^k{)up{k 2 ) 1 ^P L 1 -i- lv u a {p 2 )J':{p 1 ) {e^^-^Erfcz^ - -±-e 

\L\ L LV \L\ / y/nol 

1 \ „ „ (n\ 2 ,(j)2 



z 



e 



-iqj\L\ 



* + ^) Erfc? -^4i^ B ]})- (22) 
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where the complex vector L is defined in Eq. (|8]) and = —^5l 2 qj =F \L\/5l. Apart from 
obvious complications due to the spinorial structure of the analytic expression (|2~2"1), the 
theoretical predictions as well as conclusions concerning neutrino oscillations are very anal- 
ogous to the scalar case discussed above, and we shall not repeat these here. Instead, we 
wish to make a few clarifying remarks concerning the model under discussion. Even though 
the Gaussian-type spatial broadenings are equivalent to considering wave-packets for the 
initial (production) and final (detection) system of states, our field-theoretic approach does 
not involve explicitly wave-packets for the propagating particles. In fact, the propagation 
of the quasi-virtual intermediate states in space is consistently taken into account by their 
respective propagators without recourse to any boundary or preparation condition for the 
propagating states. Since we integrate over all times, energy conservation and flavour fac- 
torization at the source is a consequence of our model rather than an input to it. 

Because of the strict energy conservation discussed above both at production and detec- 
tion vertices, the formulas of particle oscillations in Eqs. fll3|) and fll9|) predict an infinitely 
long oscillating pattern [[II]. This is not an unexpected consequence of field theory. In 
relativistic quantum field theory, energy and three-momentum are described by commuting 
independent operators in the same way as time and space variables are. This fact is also 
manifested in the Fock-Schwinger formulation of propagators ||20|| . As a consequence, energy 
and space may in principle be measured to infinite precision, despite the spread in three- 
momentum and time for the source and target states. In fact, the solvable model under 
discussion amounts to assuming an infinite spread in time which in turn is equivalent to 
an infinite sharply peaked energy distribution. Within this framework, the use of on-shell 
kinematic conditions to predict the energy of a quantum system by its three momentum is 
not a rather well justified procedure. As we will see below, an upper limit on the length 
of particle oscillations comes from statistical (incoherent) uncertainties within the solvable 
model. 

In addition to the coherent quantum-mechanical uncertainties discussed above, one has 
to worry about statistical effects due to the energy spread of the stationary source or due 
to source's and detector's finite sizes in space which may give rise to destructive effects on 
particle oscillations. Let us discuss an illustrative example of the kind. Suppose that the 
validity conditions for oscillations in Eqs. (|14]) and (]15D are satisfied and the source being in 
a stationary situation emits neutrinos of an average energy E with dispersion AE, described 
by the function 

U(E;E,AE) = — ^— 0( -AE — \E — E\) . (23) 



AE V 2 

The function H(E) has the shape of a double-well with width AE and goes over to a delta- 
function in the limit AE — > 0. Then, for non-zero values of SE, the oscillating part of the 
averaged probability for a — > (3 transition is given by 
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/+oo 
dE U(E;E,AE)P^(E,l) 
-00 

1 r AE/2 

^ 2^e(VfV 2 ab *) — / dx cos 

V 1 AE J-AE/2 



r Am 2 / 




I- 2£ 


(-1)1 



-AE/2 

where we have implicitly assumed invariance of our model under CP, i.e., Sm(V 1 afe V 2 a6 *) = 0. 
From the last expression in Eq. ([24]), we readily see that neutrinos cease to oscillate after 
travelling a distance 

/ > L coh = L osc {E) . (25) 

nAE 

In the literature L coh is called coherence length of oscillations, which characterizes the max- 
imal scale at which neutrino oscillations can still be observed (see, e.g., ||10| ). By analogy, 
one can estimate averaging effects on oscillations due to source's and detector's finite sizes 
o x and o y , respectively. Going through a similar exercise, one finds for the oscillating part 
of the probability 

= 2 W) ^ Bin ( 2551 ) cos ( |SL ) , (26) 

where of = o~ 2 + cr 2 . Despite its different origin, the validity condition for neutrino oscilla- 
tions, first observed by Pontecorvo in [JTJ] , is then quite analogous to Eq. (0), i.e. 

o x < L osc . (27) 

As opposed to [pjj , we do not find any new upper bound on L coh within our solvable model 



if AE = and the requirements given in Eqs. (|14]) and (p7|) hold true. In fact, in such 
a case neutrinos would oscillate to infinite distances. Nevertheless, in realistic situations 
of production and detection of asymptotic states, AE is different from zero, leading to a 
physical cut-off for L coh [p],^|,p!0[16|. For example, for solar neutrinos, the energy dispersion 



AE of the production system is of the same order of the energy thermal spread AE e ~ 1 
keV of electrons captured by 7 Be. 

In summary, we have studied neutrino oscillations in space in a realistic model in which 
both the source and target states have Gaussian-form localizations. Such a model is analyt- 
ically solvable in terms of complementary error functions within the perturbative framework 
of field theory. This enables to take under close scrutiny various approximations used in the 
literature to arrive at the known formula (|l|) of neutrino oscillations, such as the steepest 
descent method or the naive pole- dominance approximation. We find that the actual space- 
dependence of particle oscillations obtained by the solvable model (c.f., Eqs. (|7|), © and 
(P2|)) is very sensitive to the localization spread 51 of the production and detection states. 
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If 81 is of microscopic size such that Sl 2 E <C I, we then recover the usual picture of particle 
oscillations through spherical waves given in Eqs. (|]) and flllf ). If, however, the distance I 
between the source and the detector is such that 5l 2 E ^> I but 51 <C L osc and I, neutrino 
oscillations occur through plane waves and their apparent oscillation length as seen by the 
detector is strongly correlated with the average three-momentum p of the initial particles 
producing neutrinos. Such a possibility may introduce ambiguities in the determination of 
the mass differences of neutrinos. Which of the two afore-mentioned forms of propagation 
dominates in different neutrino experiments and how can these be disentangled are open 
challenges that present and future experiments may have to face. 
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APPENDIX A: THREE-DIMENSIONAL INTEGRALS 



In this appendix we shall list analytic expressions of useful three-dimensional integrals oc- 
curring in the model under discussion. For this purpose, we first give the formula pertaining 
to the Gaussian integration 



d z x e~ id * e~ s ' &x = vr 3 / 2 fe 3 e-3 fa a . (A.l) 



Another useful integral is the scalar propagator in three Euclidean dimensions 

A( W(S^FWn;- (A ' 2) 

The above integral can be solved exactly by virtue of Cauchy's theorem. The intermediate 
steps of integration may be described as follows: 

\ r+oo r+1 e Aq\ lz I r+oa pMg] 1 _ g-*l?l' 

a (0 = t-? / l<?f rf l<7l / dz , ■ j- = --n 



Att 2 Jo 7-i gj — \q\ 2 + ie An 2 ! Jo q 2 — \q\ 2 + ie 

j r+oa P M?[l 

IgHgl , ,w a , .. , (A.3) 



An 2 lJ-oc q]~\q\ 2 + ie 

with \l\ — I. The integrand in the last expression of Eq. ( |A.3| ) has two complex poles situated 
at |g| = ±qj He. For the case I > 0, we may analytically continue the contour of integration 
over the upper half of the complex \q\ plane. Since the integration over the complex semi- 
circle boundary vanishes, the only non-vanishing contribution to the integral ( |A.3| ) comes 
from the pole qj + ie. Thus, we find 

A (0 = - A ( 2 ") l im (1*1 - <U - fe) 2 '^iT! • = " FT el9jl ■ (A ' 4) 
4n z l \q\->qj qj — \q\ z + iE AtiI 

In the following, we shall consider the integral contained in Eq. (||). More explicitly, we 
have 



— * 

Ai(k,p, I, Sy, Sx) 



^3- e -Sx 2 (p-q) 2 /4 e -8y 2 (k~q) 2 /4 & iql 



(2tt) 3 q] - \q\ 2 + ie 

e -¥x 2 W-\5y 2 \k\ 2 e -3( fa2+ ^ 2 ^1 2 sin(|g1|L|) 

Jo m<1] qhWTTe ' (A ' 5) 



where L is defined in Eq. (§). The integral ( |A.5|) is generally convergent since the integrand 
falls off exponentially in the limits \q\ — > ±oo. However, the theorem by Cauchy does not 
apply in this case, since the integral over the complex semi-circle with infinite radius does 
not vanish. Therefore, we must proceed differently. 

The integral we must calculate in Eq. ( |A.5|) has the following analytic form: 
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1 r+oo ryp-a^x 2 I p ixy _ p -ixy\ 

^ s «L ix — — -• (A - 6) 

where a = \^5x 2 + 5y 2 > 0, b = —iqj + e/ (2qj) (Jfteb > 0), and y = \L\ is a complex variable. 
Instead of considering the integral I(a,b,y), we shall first calculate 

1 r+oo Tp -a 2 x 2 ixy 

and then obtain I(a, b, y) by the obvious relation 

I (a, b,y) = F{a,b,y) - F{a,b,-y) . (A.8) 
We use Schwinger's representation for the integrand in Eq. (|A.7 ), viz. 

F(a,b,y) = - lim e aH2 [* * dt t / + °° dx x e'^^ e ixy . (A.9) 

l R-*+oo J a Jo 

In order to avoid possible singularities originating from the complex pole (x 2 + 6 2 ) -1 of the 
propagator, we could perform the usual Wick rotation and regard x as a pure imaginary 
variable. We shall no longer dwell upon this point, since it also occurs even if a — > 0, in 
which case the result may be read off from Eq. ( |A.4j ). Thus, we can formally carry out the 
Gaussian integration in (|A.9|) and rewrite F(a, b, y) in the following forms: 

F(a, b,y) = -- lim e aH2 dt t / + °° dx e^ 2 e lxy 

2 R-^+oo J a dy J-oo 

lim e a2b2 f R/b d lye-^ 2 -y 2 '^ 

R^+oo J a t 2 

lim e a2b2 e-* % y^-vl^ . (A.10) 

R-+oo Ja t 2 y ^ ' 



Making the variable substitution z = tb — y/{2t) and taking the limit R — > +oo where 
appropriate, we finally arrive at 

F(a,b,y) = - ^ 2 Erfcfafe- 1-) _ b -^l l im dt e -^-vV(^) 

v y> 4 V 2a' 4 R-+oo J a 

(A.ll) 

In Eq. ( |A.11| ), Erfc(z) is the complementary error function of a complex argument defined 



as 



Erfc(z) = -= / dt e~ t2 = 1 = / dt e~ t2 . (A. 12 

V7T Jz \l 71 Jo 



The last /^-dependent term on the RHS of Eq. (|A.11| ) is symmetric under the interchange 
y — > —y and cancels in I(a,b,y) when the difference in Eq. (|A.8|) is formed. The so-derived 
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analytic result agrees with that of Ref. for y > which has been extended here to include 
complex values of y. After restoring all the terms present in Eq. (|A.5| ), we eventually obtain 



A4(k,p, 1, 5y, 5x) 



8tt\L\ 



e _I fa a ( |ji|a + ^ ) _l^ ( | fc |a + ,a ) 



X 



e iii\L\ Erfc(-|% - 13 



(A.13) 



with 5l 2 = 5x 2 + by 1 . 

It is useful to list a few values of Erfc(z) for special z's as well as derive appropriate 
expansions of its argument which have been used extensively throughout our discussion. In 
particular, we have 



Erfc(+oo) = 0, Erfc(0) = 1, Erfc(-oo) 



(A. 14) 



The asymptotic behaviour of Erfc(2;) for \z\ 1 with $lz > [2~2"| may be obtained by a 
series expansion, viz. 



Erfc(z) 



1 x 3 x ... (2m - 1) 



m=l 



(2z 



2\m 



If ^tz < 0, one may employ the following property: 

Erfc(z) = 2 - Erfc(-z) , 



(A.15) 



(A.16) 



and then apply the asymptotic expansion in Eq. ( |A.15| ). By means of Eqs. ( |A.15| ) and 
( |A.16| ), we find that in the limit 51 — > 0, the expression in Eq. (|A.13|) becomes identical 
to that given in Eq. ( |A.4|) . Furthermore, if 5l 2 qj , 5l 2 \k\ <ti I with k = p, we obtain the 
asymptotic expression in Eq. (|I2|), which is in agreement with results derived in ]11[ from 
different considerations. The latter constitutes another non-trivial check of the validity 
of our analytic results. On the other hand, if 5l 2 qj , 5l 2 \k\ ^> I but 51 <C / , L osc , the 
transition amplitude M. has a non-trivial dependence on the Erfc functions in Eq. ( |A.13| ) 
which multiply the exponential factors exp(iqj\L\) and exp(— iqj\L\). However, the term 
depending on exp(iqj\L\) leads to a Gaussian-like function exp[— <5/ 2 (|A;|— qj) 2 /4] that controls 
momentum conservation, whereas the term proportional to exp(— iqj\L\) turns out to be 
extremely suppressed by a factor exp[— {1/51) 2 }. As a consequence, we find that up to an 
overall normalization constant Af, the amplitude M. exhibits a plane-wave behaviour in this 
kinematic range (see also Eq. flT7|)). 

Finally, it is amusing to discuss the S-matrix limit of formula ( [A.13j ). As can be seen 
from Eq. (0), this corresponds to the case 5x, 5y I. For simplicity, we work in the k = p 
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approximation and keep only the leading term in the asymptotic expansion (|A.15|) . After 
some algebra, we find for kl > 



M{k = p,l,5l) 



1 



4tt\L\ 
1 

5/ 3 7T 3 / 2 



e -\SP(\kf+q]) e i qj \L\ + 



51 e~ l2 / &l2 e m 



4tt 3 / 2 p + Sl*(q] - \k\ 2 )/4 - iSPkl 



i7ie i qj W\k\ si 



\k\ 2^? 



e -\8l*{\k\- q] f 



ikl 



+ v- 



+ me ikl 5+{q 2 - \k\' 2 ) 



q 2 - \k\ 2 

where V stands for principle value. Recognizing now the fact that 



(A.17) 



S(x) = lim 4= e"" 2 " 2 



(A.18) 



and noticing the cancellation of the additional delta function S + (q 2 — \k\ 2 ) in Eq. (|A.17 ), it 
is then not difficult to reproduce the known S-matrix result 



M{k = p, I, SI — >■ oo) 



Akl 



5PtT 3 / 2 



in5+(q 2 - \k\ 2 ) + V 



(A.19) 



In the usual S-matrix treatment, one first considers the limit 8x, 5y — > oo in Eq. (|6|) and 
then integrates over the corresponding three-momentum delta functions. Up to an irrelevant 
normalization, the result thus obtained is identical to that of Eq. ( |A.19[) . 
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